Optimized online filter stack spectrometer for ultrashort X-ray pulses* 
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Currently, with the advent of high-repetition-rate laser-plasma experiments, the demand for online diagnosis 
for the X-ray spectrum is increasing because the laser-plasma-generated X-ray spectrum is very important for 
characterizing electron dynamics and applications. In this study, scintillators and silicon PIN (P-type-Intrinsic- 
N-type semiconductor) diodes were used to construct a wideband online filter stack spectrometer. The X-ray 
sensor and filter arrangement was optimized using a genetic algorithm to minimize the condition number of the 
response matrix. Consequently, the unfolding error was significantly reduced based on numerical experiments. 
The detector responses were quantitatively calibrated by irradiating the scintillator and PIN diode with various 
nuclides and comparing the measured y-ray peaks. A prototype 15-channel spectrometer was developed by 
integrating an X-ray detector with front- and back-end electronics. The prototype spectrometer could record X- 
ray pulse signals at a repetition rate of 1 kHz. Furthermore, an optimized spectrometer was employed to record 
the real-time spectra of laser-driven bremsstrahlung sources. This optimized spectrometer offers a compact 
solution for spectrum diagnostics of ultrashort X-ray pulses, exhibiting improved accuracy in terms of spectrum 
measurements and repetition rates, and could be widely used in next-generation high-repetition-rate high-power 


laser facilities. 
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I. INTRODUCTION 


Laser-driven plasma-based electron accelerators and X- 
ray sources powered by ultra-intense laser technology have 
been extensively researched in recent years. With an ac- 
celerating gradient of approximately 100 GeV/m in laser 
plasma accelerators, electrons can be accelerated to hundreds 
of MeV within a few millimeters [1]. Subsequently, these 
electrons can generate X-rays through betatron radiation [2], 
inverse Compton scattering [3-5], and bremsstrahlung [6], 
thus providing tabletop alternatives to large-scale conven- 
tional accelerator-based X-ray sources. These X-ray sources 
have advantages such as femtosecond duration, micron-sized 
source dimensions, and a wide spectral range [7], thereby ex- 
hibiting tremendous potential for multiple applications [8], 
including biological radiography [9], non-destructive testing 
[10, 11], and high-energy-density physics [12, 13]. 

The research and application of laser-driven tabletop X-ray 
sources require a unique set of diagnostics [14]. Among them, 
X-ray spectrometers are particularly important as they char- 
acterize electron dynamics in plasma [4, 5, 15, 16] and are 
useful in various applications [13]. Multiple types of spec- 
trometers have been developed to cover a wide range of pho- 
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ton energies, from a few keV to tens of MeV, in laser-plasma 
experiments. For X-rays below 30 keV, charge-coupled de- 
vices (CCDs) operating under single-photon counting mode 
provide high-resolution X-ray spectra by identifying individ- 
ual X-ray photon signals in a single image frame [17, 18]. 
Crystal spectrometers can be employed for X-rays below 100 
keV, and the spectrum can be measured by detecting the an- 
gular distribution of X-rays diffracted by the crystal [19, 20]. 
For X-rays with energies in the range of hundreds of keV, 
an electron-track-based spectrometer has been proposed that 
reconstructs X-ray energies by detecting Compton electron 
tracks and energies using silicon trackers [21]. For X-rays 
greater than 1 MeV, the Compton spectrometer converts the 
X-rays into electrons via forward Compton scattering. Comp- 
ton electron spectra are measured using a magnetic spectrom- 
eter, from which X-ray spectra can be derived [22-26]. Al- 
though the aforementioned spectrometers are sensitive to X- 
rays within relatively narrow energy ranges, the filter stack 
spectrometer (FSS) enables the diagnosis of X-ray spectra 
over a wide band, typically ranging from tens of keV to hun- 
dreds of MeV. A FSS employs a stack of filters interlaced 
with X-ray sensors [27—29], and the characteristic parame- 
ters of the X-ray spectra can be reconstructed from the re- 
sponse matrix (RM) and sensor signals. Additionally, FSSs 
have the advantages of simple manufacturing and simple op- 
eration; therefore, they have been widely used in laser-plasma 
experiments. 


Currently, many FSSs utilize image plates (IPs) as X-ray 
sensors [27-29] because of their reusability, high sensitivity 
to X-rays, versatility, and resistance to electromagnetic pulses 
[30, 31]. However, the IP signals must be read offline using an 
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IP scanner, which takes approximately 10 min for data read- 
ing. This offline readout time does not meet the requirements 
of high-repetition-rate laser-plasma experiments. Given the 
importance of high repetition rates in studying experimental 
laws and optimizing accelerators and X-ray sources [32-34], 
there is a need for online diagnosis of the wideband spectrum 
of laser-induced X-rays. 


Scintillators and semiconductors are commonly used in 
X-ray sensors [35]. Scintillators convert X-rays into visi- 
ble scintillation light, which can be easily transformed into 
electric signals using a photoelectric converter. In contrast, 
semiconductors convert X-rays into charge carriers (electrons 
and holes). Electric signals are generated at the electrode 
as the charge carriers drift under a bias voltage. The de- 
cay time of scintillators and the charge carrier drift time of 
semiconductors typically range from a few to hundreds of 
nanoseconds. Therefore, electrical signals can be recorded 
and processed online using an electric readout system. Con- 
sequently, an FSS constructed with scintillators or semicon- 
ductors would allow for a fast repetition rate, aligned with 
the current state-of-the-art and next-generation laser systems. 
Recently, online FSSs based on scintillators have been devel- 
oped [36-39] utilizing CCDs or complementary metal-oxide- 
semiconductor (CMOS) cameras as photoelectric converters. 
Notably, scintillator-based FSSs have already been employed 
in experiments [40-42], demonstrating exceptional perfor- 
mance in wideband X/y-ray detection. However, certain as- 
pects of the currently developed online FSS techniques re- 
quire further improvement. 


First, the light collection efficiency is low because the lens 
of CCDs or CMOS cameras covers only a small solid angle. 
This limitation can result in high statistical uncertainties, par- 
ticularly under conditions of low X-ray fluence. Second, the 
calibration and stability of the scintillation-light-collection ef- 
ficiency pose challenges owing to the non-solid nature of light 
couplings. Finally, the FSS configuration has not been ade- 
quately optimized, resulting in a severely ill-conditioned RM. 
These limitations may contribute to a relatively large error in 
X-ray spectrum determination. 


In this paper, we propose an optimized online fast X-ray 
spectroscopy system (FSS) and present its implementation 
herein. To address the ill-conditioned RM, we employed a 
genetic algorithm [43] to optimize the sensor and filter ar- 
rangement. To achieve robust light coupling and a high light 
collection efficiency, we utilized PIN diodes with large sensi- 
tive areas and low noise as the photoelectric converters. Fur- 
thermore, PIN diodes were employed as sensors in the chan- 
nels closest to the X-ray sources, enabling the direct detection 
of low-energy X-rays and improving the unfolding accuracy 
within an energy range of tens to one hundred keV. The pro- 
posed configuration allows for direct calibration of the light- 
collection efficiency and alleviates the ill-conditioned nature 
of the RM. According to numerical experiments, the relative 
deviation of the true and unfolded spectra can be decreased 
to approximately 16%. The electrical signals from the PIN 
diodes are amplified using in-house built trans-impedance 
amplifiers (TIAs) and main amplifiers, and then digitized by 
a multi-channel data acquisition (DAQ) system based on the 


112 
113 
114 
115 
116 
117 
118 
119 
120 


12 


122 


123 


124 


125 


126 
127 
128 
129 
130 


13 


132 
133 
134 
135 
136 
137 
138 
139 
140 


14 


142 


143 
144 


145 
146 
147 
148 
149 


150 


151 


152 
153 
154 
155 
156 
157 
158 
159 
160 


16 


DRS4 chip [44]. The proposed online FSS integrates filters, 
X-ray sensors, front-end electronics (FEE), and the DAQ sys- 
tem, resulting in a compact and user-friendly system. No- 
tably, the proposed online FSS exhibits sensitivity to X-rays 
ranging from tens of keV to hundreds of MeV, enabling the 
recording of X-ray pulse signals at a repetition rate of 1 kHz. 
By contrast, CCD- or CMOS-based FSSs typically have a 
maximum repetition rate of 100 Hz [36]. Furthermore, the 
proposed online FSS demonstrates an energy measurement 
error of less than 1%, whereas traditional IP systems have 
a 20% energy measurement error [29], thereby providing a 
more accurate unfolding spectrum. 


I. METHODOLOGY 


A. Spectrometer optimization 


The most significant drawback of the FSS is the ill-posed 
linear system, which results in substantial errors in spectrum 
measurements. Because the scintillator can function as both 
filter and sensor, parameters such as the layer thickness, num- 
ber of layers, filter density, and scintillator density need to 
be finely tuned to alleviate the ill-conditioned nature of the 
RM. Rusby et al. [36] optimized their online FSS based 
on the difference in the scintillator output within the energy 
range of interest. The optimal scintillator should exhibit the 
largest difference in scintillator output. However, this opti- 
mization criterion was spectrum-dependent and the scintilla- 
tor thicknesses were kept constant, limiting the optimization 
efficiency. 

To propose a universal optimization method, we first es- 
tablish a formal model for encoding an X-ray spectrum S (E) 
into the experimental channel D; for an N-channel FSS. The 
mathematical definitions and notation employed in this study 
are consistent with those used by Fehl et al. [45, 46]. D; can 
be written as follows: 


Emax 
0 
(1) 


where D; represents the experimental channel data, 
R;,(E)i = 1% denotes the response function, di represents 
the noise-free channel data, and e€; represents uncertainties 
and noise. Equation (1) can be expressed in vector-matrix 
notation as follows: 


D=RS+AD=d+AD. (2) 

The errors that affect the unfolding accuracy can be di- 
vided into two parts. The first part is the perturbation AD 
superimposed on D during the data-gathering process, which 
includes factors such as statistical uncertainty, signal noise, 
digitization errors, and uncorrected signal baselines, among 
others. The second part is the bias between the experimen- 
tally measured or simulated RM and the actual RM, expressed 
as R — R + AR. The presence of AR can be attributed to 
drifts and uncertainties in the fitting parameters, such as the 
calibration of the light yield and light collection efficiency 


12 discussed in Sect. IC. Because (AR)S is analogous to AD 
1e3 [46], only AD is discussed here. 

Solving Eq. (2) poses challenges because it is an ill-posed 
165 problem. Various unfolding algorithms can be employed to 
address this problem, including regularization methods, least- 
squares spectrum adjustment, parameter estimation, iterative 
unfolding methods, and the maximum entropy principle [47]. 
However, these unfolding algorithms often struggle to pro- 
vide mathematically rigorous or realistic estimates of the er- 
171 ror propagation relationship between the spectral unfolding 
error AS and the data perturbation AD. Nevertheless, it is 
worth noting that the error propagation relationship between 
AS and AD is primarily determined via the response func- 
tion R. For instance, an upper-bound measure for AS can be 
estimated from AD and the condition number of the RM is 
denoted by cond(R), as described in [45]. 


||AS||/||S|| < cond(R)|[AD]|/||D]], 
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(3) 


where ||D|| = can D?)2 represents the norm of D. RMR 
can be either a square or non-square matrix [59], and the con- 
dition number cond(R) can be calculated via singular value 
decomposition. By reducing cond(R), for example, by al- 
leviating the ill-conditioned nature of the RM, the upper- 
184 bound measure of AS can be lowered, resulting in improved 
precision in terms of the spectrum measurements. Further- 
more, cond(R) is determined solely by the FSS configuration. 
Hence, cond(R) can serve as a measure of merit for FSS op- 
timization, independent of the spectrum. 

Another challenge in spectrometer optimization is the 
multi-parameter optimization problem. FSS typically com- 
prises numerous filter and sensor layers, and the thickness 
and density (determined by the material) of each layer can 
affect cond(R). Consequently, multiple parameters must be 
optimized. Owing to the practical limitations on the layer 
thickness for both the filter and sensor, as well as the absence 
of derivatives and linearity in this problem, the genetic algo- 
197 rithm is a suitable method for optimization. 

Several factors restrict the range of the optimization pa- 
rameters. In an FSS, filters and sensors closer to the X-ray 
source are typically thinner, whereas those farther away are 
thicker. This design ensures that most low-energy X-rays 
deposit their energy in the layers near the sources and high- 
energy X-rays do so in the layers farther away. Consequently, 
RM R has a smaller condition number (cond(R)). To address 
the low light-collection efficiency of thin scintillators, which 
can result in significant statistical errors, sensors with a 300 
jum depletion layer thickness, such as PIN diodes, were used 
in the layers near the X-ray source. For custom production 
convenience, the scintillators are made of Gd3AloGa30 9 
(GAGG) [48], while the filters are made of aluminum or cop- 
211 per. After determining the number of layers, cond(R) should 
212 be minimized because it is scale-dependent. Considering that 
the number of layers is a tradeoff between the applicable en- 
ergy range and the complexity of the signal readout system, 
15 layers were selected, including six PIN diodes and nine 
scintillators. The overall configuration of the proposed online 
217 FSS is listed in Table 1. The layer number represents the dis- 
218 tance from the X-ray source with the filter closer to the source 


185 


in each layer. As an aluminum film of 50 um is always used 
as the first layer for electromagnetic shielding, the parameters 
to be optimized are the filter thickness in layers 2—15 and the 
scintillator thickness in layers 7-15. 


Table 1. General configuration of the proposed online FSS 


Layer Filter Sensor Layer Filter Sensor 


1 Al Si-PIN 9 Cu GAGG 
2 Al Si-PIN 10 Cu GAGG 
3 Al Si-PIN 11 Cu GAGG 
4 Al Si-PIN 12 Cu GAGG 
5 Al Si-PIN 13 Cu GAGG 
6 Al Si-PIN 14 Cu GAGG 
7 Cu GAGG 15 Cu GAGG 
8 Cu GAGG 


A key aspect for implementing the genetic algorithm is 
the fitness assignment method. As the population consists 
of 2000 individuals and there are 300 iterations, approxi- 
mately 6 x 10° cond(R) calculations, i.e., RM calculations, 
are required. Accurately RM calculations using the Monte 
Carlo method [28, 36, 37] consume significant computational 
229 power. To reduce the computational power requirement, a 
simplified RM calculation model that utilizes the X-ray mass 
energy absorption coefficients and mass attenuation coeffi- 
cient is proposed. First, the expected energy deposition of 
a single X-ray with energy ŒE in the i-th sensor, denoted by 
Edep,i(E), can be calculated as 


i 


i-1 
Edep (E) = (1—e~Hen (Esi). TT] eti (E)s; JI eH (ENG E, 


j=l 

(4) 

where Hen i (E) represents the X-ray mass energy absorption 

coefficient of the i-th sensor, p; (E) represents the X-ray mass 

attenuation coefficient of the j-th sensor or filter, sj and fj 

denote the thickness of the 7-th sensor and filter, respectively. 
Thus, the RM equation can be written as 


j=1 


R,(E) = Eaep,i(E) - ESC). (5) 


where ESC; are the energy-signal coefficients of the i-th 
channel. For the PIN diode, 


where 3.62 eV is the pair creation energy of silicon [49], q is 
the elementary charge, and p is the gain of the electronics in 
ohm. For GAGG, 


ESC; = (LY x LCE x QE); x q X pi, (7) 


where LY is the light yield of GAGG, LCE is the light col- 
lection efficiency, and QF is the quantum efficiency of the 
photon detector. If we disregard the nonlinearity of the scin- 
tillator light yield and electronics, ESC; becomes indepen- 
253 dent of th X-ray energy FE. Therefore, the practical response 
254 equation to be solved can be derived using Eqs. (1) and (5) as 
255 follows: 


D/ESC = 


256 dep S unfold: (8) 


257 Where Sunfola is the spectral unfolding and Eqep is the matrix 
notation of Egep,;(£). An error AS exists between Sunfola 
and the real X-ray spectrum S owing to perturbation AD in 
D. In Eq. (8), Egep is equivalent to RM R, and can be deter- 
mined from Eq. (4). Therefore, the optimization process for 
reducing cond(Ege,) is illustrated in Fig. 1. 

The genetic algorithm was implemented using the Sheffield 
Genetic Algorithm toolbox [50]. The population consisted 
of 2000 individuals, with each individual containing 14 fil- 
ter thicknesses and nine scintillator thicknesses. These thick- 
nesses are represented as floating-point numbers and encoded 
as bit strings to create chromosomes. Each individual con- 
tained 24 binary chromosomes, representing the 14 filter 
thicknesses and nine scintillator thicknesses. To consider 
practical layer thickness limitations, lower and upper limits 
were set for the filters and scintillators, 0.01 and 50 mm and 
2 and 50 mm, respectively. Each binary chromosome uses 20 
bits to ensure sufficient precision. For the iterations, all filter 
thicknesses were set to zero, and any thickness below 0.1 mm 
was also set to zero. The generation gap, mutation probabil- 
ity, and crossover probability were set to 0.9, 0.01, and 0.7, 
respectively. The termination criterion was reached when the 
iteration time was 300. 

The minimum and average Cond(Egep) values during the 
iterations are shown in Fig. 2. The optimal Cond(Egep) is 
approximately 4 x 104, which represents a significant reduc- 
283 tion compared with the unoptimized online FSS configura- 
234 tion. For instance, the previously reported CsI FSS [37] ex- 
285 hibited a Cond(Eqep) of 2 x 101? (considering 15 channels). 
236 The optimal configuration obtained during the iterations is 
27 listed in Table 2. The 0.01-mm precision is maintained, con- 
288 sidering that the machining tolerance of the GAGG crystal is 
ze9 0.02 mm. 
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Table 2. Optimized configuration of the proposed online FSS 


Layer Filter Thickness Sensor Thickness 
(mm) (mm) 
1 Al 0.05 Si-PIN 0.3 
2 Al 0 Si-PIN 0.3 
3 Al 0 Si-PIN 0.3 
4 Al 2.82 Si-PIN 0.3 
5 Al 9.95 Si-PIN 0.3 
6 Al 24.34 Si-PIN 0.3 
7 Cu 3.53 GAGG 2 
8 Cu 10.07 GAGG 2 
9 Cu 4.73 GAGG 4.44 
10 Cu 6.97 GAGG 13.04 
11 Cu 0.57 GAGG 12.15 
12 Cu 13 GAGG 33.74 
13 Cu 1.14 GAGG 21.65 
14 Cu 18.24 GAGG 33.2 
15 Cu 34.52 GAGG 50 


20 Itis worth noting that the optimization process using a sim- 
21 ple RM calculation model has certain limitations. First, be- 
2932 cause the mass attenuation and mass X-ray absorption coeffi- 
23 cients are available only within an energy range of 1 keV to 20 
204 MeV, the optimization process considers X-rays ranging from 
25 10 keV to 10 MeV. Second, it is well known that the simple 


26 RM calculation model does not consider secondary effects, 
297 Such as Compton scattering and pair production processes, 
which become dominant in photon-matter interactions start- 
ing at hundreds of keV. Consequently, the simple RM calcu- 
lation model results deviate from the actual values above hun- 
dreds of keV. However, it still provides a reasonable estima- 
tion of the actual RM, and we demonstrated that optimization 
is an effective approach for achieving a partially optimized 
structure. This is evident from the significant reduction in the 
unfolding error, as demonstrated next. 


298 
299 
300 


30 


302 
303 
304 


305 


306 B. Unfolding spectrometer response 


Numerical experiments were conducted to test the unfold- 
ing procedure and evaluate the accuracy of the unfolded radi- 
ation spectra. Using Eqs. (2) and (5), the expected energy de- 
position Egep, X-ray spectrum S, energy-signal coefficients 
ESC, and perturbation AD were calculated and modeled to 
simulate the channel data. 

First, the optimized online FSS was modeled using 
GEANT4 [51], a simulation platform using the Monte Carlo 
method, to obtain an accurate Eg.,. The simulation model 
consists of the filters, sensors, scintillation photon detector, 
PIN diode carrier board, collimator, shielding, and other me- 
chanical parts. The incident X-rays were simulated as pencil- 
like beams uniformly and randomly distributed within a 6- 
mm_-diameter circle, matching the collimator diameter. The 
energy range of the incident X-rays was set from 10 keV 
to 200 MeV. For each energy bin, the X-ray energies were 
uniformly and randomly sampled and 10° X-rays were simu- 
lated. The energy deposition curves of the X-ray beams im- 
ses pinging on the online FSS, i.e., Egep, are shown in Fig. 3. 
The X-ray spectrum S is modeled as follows: 


S(E) = SBetal E) + SBrems(£) + Spa(£), 
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(9) 


where Sgeta(E) = AE°exp(—(E) represents the on-axis 
329 betatron component, Sprems = Biv — In(£)) represents 
the on-axis bremsstrahlung component, and Sga(E) 
Cexp(—1E) represents the background component [28, 52]. 
Note that E is in MeV. Parameters A, a, 6, B, v, C, and 7 
define the shape of the curves and the fluence ratios of the 
three components. Specifically, they are set to 12.47, 0.768, 
7.52, 4.09 x 1074, 3.88, 0.01, and 1, respectively. In this 
scenario, the betatron radiation critical and peak energies are 
0.25 and 0.1 MeV, respectively, and the fluence ratio of the 
background component is 7%. This represents a typical be- 
tatron radiation spectrum with background [28]. Assuming 
X-ray energies ranging from 10 keV to 10 MeV and a to- 
tal photon count of 10°, the simulated energy depositions in 
the online FSS channels were obtained by multiplying the ex- 
pected energy deposition Ege, by the photon count in each 
energy bin, S(E)dE. 

34 ESC can be estimated using Eqs.(6) and (7). GAGG-HL- 
34 type GAGG crystals manufactured by EPIC CRYSTAL Co., 
347 Ltd., with a light yield LY of 54000 photons/MeV were used 
34s in this study. The light-collection efficiency LCE was ap- 
3 proximately 60% under typical packaging and optical cou- 
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Population Initialization Fitness Assignment Termination 
Generate a number of random FSS Calculate E4ep and Cond(Egep) for ee 
configurations in the given ranges each configurations criteria met? 


Return the FSS 
configuration with 
minimum Cond(Egep) in 


Generate all generations 
New Crossover 
Population 


Fig. 1. Optimization process of the proposed online FSS using the genetic algorithm 


10" OMinimum cond\E,,,) 32 efficiency QE of PIN diodes, such as the 0.3-mm Si-PIN, 
10° Ex is 33 is approximately 90% at a wavelength of 530 nm. Because 

E Average cond(E sep) f 3 R 
io EA 34 ESC may vary owing to the different scintillators, packag- 
~ FA 355 ing, and optical couplings, it was experimentally calibrated 
rye 10 Lx 356 as described in Sect. II C. The FEE gains (p) are the product 

5 E À 357 Of the gains of the TIAs and main amplifiers, which collec- 

S 10 E04 ; ; ; 

Oo  Eo4 358 tively determine the overall gain of the electronic unit. The 
10° 359 gains of the TIAs and main amplifiers were determined indi- 
ie 360 Vidually to ensure that the signal amplitudes conformed to the 

E 361 dynamic ranges of each stage. A detailed discussion of p is 
| | 


o 50.100 750 “500° “250. “300. see provided in Sect. NIB and the values are listed in Table 4. 

generators 33 The perturbation AD is primarily caused by the statistical 
uncertainty o (E) s+, where E represents the energy deposition 
35 in the detector, and the electronic noise Cnoise. The total per- 
3 turbation Gtotal is calculated as Ctotal = o(E)?, +02 


noise’ 


37 The statistical uncertainty in the PIN diode channels can be 


36: 


P 


Fig. 2. (Color online) Evolution of Cond (E4ep) during the genetic 
algorithm iterations 


Channels GAGA A aes estimated using o(E)st = 2.355VF - E-W, where F = 
crena dace |E w 369 0.12 is the Fano factor and W = 3.62 eV is the energy re- 
e a 2 370 quired for the formation of a charge carrier pair. Statistical 
Channel10 GAGG 1o25 371 uncertainties in the GAGG channels were estimated based on 
ines 3 372 a previous calibration experiment conducted in our labora- 
Channel? GAGG 1038 373 tory using radioactive sources [53]. In particular, o(E)st = 
Channelé Si-PIN a 3 . Z : 

Channels Si-PIN 5 az 0.09 - 662 x 10° - 4/ 662x107 The electronic noise model is 
N 1043 37s complex and beyond the scope of this study. The calculation 
anaes 8 376 Method for the electronic noise in a high-speed waveform- 

" 10°5 


1 


i ið 10? 377 sampling detector is described in detail in [54]. According to 

Energy (MeV) sve Cang et al. [54], when the pulse amplitude fully utilizes the 

39 ADC range and the integration time window is 150 ns, the 

Fig. 3. (Color online) Response matrix expressed as Eacp obtained 380 electronic noise reaches approximately 0.12%. Considering 

via Monte Carlo simulations 381 that our time integration window is 600 ns and the pulse am- 

332 plitude may vary from pulse to pulse, a 1% electronic noise 

333 error estimation is reasonable. Therefore, an electronic noise 

350 pling conditions [53]. The maximum emission wavelength 384 Of Cnoise = 0.01 x D; is applied to all channels, where D; 
351 for GAGG is approximately 530 nm [48], and the quantum sss represents the noise-free channel data. 


Based on the simulated channel data, Dg;,,, the unfolded 
3837 X-ray spectrum was obtained by solving Eq. (8) using the 
388 EXpectation-maximization method [55]. The convergence de- 
389 gree is indicated by the normalized mean absolute distance 
300 (MAD) between the fitting and simulated channel data. 


386 


391 (10) 


1 |D; fit — Di simi 
MADaata = : EAN 
dat N D Di sim 


32 where N denotes the number of channels. Using a constant 
33 function of unity as the initial guess for Sunfoa(E)dE, the 
324 fitting errors for different numbers of iterations are plotted in 
Fig. 4(a). In the case of the optimized GAGG array (repre- 
3 sented by the red open circles), the fitting error initially de- 
397 creases rapidly to a low level within the first few steps and 
e eventually reaches < 1% after 50 steps. The deviation be- 
tween the unfolded spectrum and the spectrum model func- 
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4v tion is also indicated by the normalized MAD between them. 
1 |[Sunfoa (E)dE — S(E\dE| 
MAD sec = marae 
p M >. S(E)dE ( ) 


42 where M denotes the number of energy bins. As shown in 
43 Fig. 4(b), the divergence can reach a minimum of approxi- 
44 mately 16% in an energy range of 10 keV to 10 MeV when 
s the fitting converges. In this study, 100 iteration steps were 
chosen, which is sufficiently large to observe the convergence 
47 properties, obtain a reasonable spectrum shape, and fit the re- 
sponse of the FSS, as shown in Fig. 5. It should be noted 
that the choice of the unfolding algorithm can significantly af- 
fect the unfolding results, and the expectation-maximization 
method was selected because of its ability to avoid the need 
for many additional constraints and its greater universality. 
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Fig. 4. (Color online) (a) Converge curve and (b) unfolding error for 
three kinds of online FSSs 


“3 For a comparative analysis, the CsI array proposed by 
44 Behm et al. [37] and the BGO array proposed by Rusby et 
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Fig. 5. (Color online) Fitting data (right) plotted with the corre- 
sponding spectra (left) 


al. [36] were also simulated following the aforementioned 
procedure. The channel number and noise level were kept 
the same as those in the GAGG array, that is, 15 channels. 
The thickness of the BGO crystals was 0.2 cm. The condi- 
tion numbers for the CsI and BGO arrays were 2 x 101? and 
1 x 10/2, respectively. The corresponding minimum fitting 
errors were 115% for the CsI array and 133% for the BGO 
array, respectively, as shown in Fig. 4. The unfolding er- 
ror of the optimized GAGG array was significantly reduced 
compared to that of the CsI and BGO arrays with uniform 
thickness. 


C. Quantitative calibration of the spectrometer response 


An accurate RM determination is of vital importance for 
precise spectrum unfolding. In Eq. (5), Edep,i(£) can be pre- 
cisely obtained through Monte Carlo simulations [28, 36, 37]. 
However, ESC must be experimentally calibrated for scintil- 
lator channels because the light yield LY and light collection 
efficiency LCE can vary depending on the scintillators, pack- 
ages, and optical couplings, making it challenging to model 
and calculate them. 

The most challenging aspect of the experimental calibra- 
tion campaign for LY and LCE is determining the energy 
deposition in the scintillators. Behm et al. [37] performed an 
experimental calibration by measuring the scintillation sig- 
nal resulting from bremsstrahlung interactions. They theoret- 
ically calculated the bremsstrahlung X-ray beam in GEANT4 
by simulating the collision of a typical electron beam with a 
9-mm-thick piece of lead. The energy deposition in the CsI 
array was then determined using GEANT4. Rusby et al. [36] 
calibrated their online FSS by exposing the detector to radi- 
ation sources for an extended period of time and integrating 
the camera images. Energy deposition during this period was 
also determined via GEANT4 simulation. 

In the optimized online FSS, a PIN diode with a 0.3-mm 
Si-PIN was also used as the scintillation photon detector. The 
300-um depletion layer of the 0.3-mm Si-PIN allows us to ob- 
tain LY x LCE x QE by comparing the full-energy peaks of 
the radioactive sources. One source directly irradiates the PIN 
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diode 0.3-mm Si-PIN, while the other irradiates the GAGG 
scintillator [56]. The energy deposition in the scintillators, 
that is, the gamma ray lines of the radioactive sources, can be 
accurately determined, thereby providing more reliable cali- 
bration results. 

The calibration experiment utilized 7+! Am with a 59.5-keV 
line and ??Na with a 1274.5-keV line. Charge-sensitive pre- 
amplifiers were used instead of TIAs to amplify the signals 
induced by the radioactive sources. An example of the spec- 
trum measured during the experiment for a single channel is 
shown in Fig. 6. The peak positions were obtained by fitting 
the experimental spectrum using a Gaussian curve superim- 
posed on a second-order polynomial background, denoted by 
the green lines; the background is represented by the brown 
lines. 
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Fig. 6. (Color online) Spectra obtained during the calibration exper- 
iment for one channel 


As the readout electronics are identical for GAGG and PIN 
diodes, according to Eqs. (5)-(7), LY x LCE x QE can be 
calculated as follows: 


Poaaac Edep,PIN 


LY x LCE x QE = , 
Q 3.62[eV] - Pern Edep,Gaaa 


(12) 


where Paaca and Ppr are the fitting peak positions in the 
experimental spectrum, Egep,gaca and Eqep,pIn are the en- 
ergy depositions in GAGG and PIN diodes, which are 1274.5 
and 59.5 keV, respectively. The corresponding spectra are 
shown in Fig. 6, Paaag and Pprn are 1.0646 + 0.0002 and 
0.4314+0.0001 V, respectively; therefore, LY x LCE x QE 
= 31826 + 5 electrons/MeV. 


Tl. IMPLEMENTATION 


An online FSS prototype was built and calibrated using the 
aforementioned methods. A schematic of the online FSS pro- 
totype is shown in Fig. 7. Similar to classic radiation detec- 
tors [57], the FSS prototype comprises sensors, FEE, a DAQ 
system, and various software. The prototype comprises six 
PIN diode channels and nine GAGG scintillator channels as 
X-ray sensors, which are interleaved with aluminum and cop- 
per filters. The collimator and laser sight are located in front 
of the PIN diode sensors, and the laser sight is replaced with 
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an outer collimator once the FSS is mounted. The scintilla- 
tors are coupled with PIN diodes as photon detectors. The 
signal output pins of the PIN diodes are connected to an FEE 
board. The following DAQ board provides the capability for 
16-channel waveform digitization and communication. In ad- 
dition, a power board supplies low-noise +5 V and 70 V 
power to the FEE and PIN diodes. To avoid interference from 
scattered X-rays, the filters and X-ray sensors are shielded by 
a 1-cm-thick lead layer and a 1-cm-diameter aperture allow- 
ing for X-ray incidence. The aperture is covered with a 50- 
„m-thick aluminum foil for electromagnetic shielding. The 
electronics are housed in an enclosed copper and aluminum 
box for electromagnetic shielding. The external interface in- 
cludes only a 12-V power supply cable and a USB communi- 
cation cable. 
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Fig. 7. (Color online) Cross-sectional view of the online FSS proto- 
type 


A. Filter and sensor unit 


GAGG scintillators were chosen because of their good 
mechanical characteristics, non-hygroscopic properties, high 
light yield (30-70 ph/keV), and fast decay times (~100 ns). 
The scintillators had front faces measuring 1 cm x 1 cm 
with variable thicknesses. The filter and scintillator config- 
urations are listed in Table 3. The prototype configuration 
was obtained using the original version of our optimization 
code without proper constraints on the thickness, resulting in 
a configuration different from that listed in Table 2. Neverthe- 
less, the condition number of the prototype response matrix 
is 6 x 10°, which is close to the optimized condition number 
presented in Sect. II A (compared with the uniform-thickness 
configuration), and an improved performance is expected. 

The reflection layers of the scintillators are enhanced spec- 
ular reflectors (ESRs), 65-um polymers with high reflectance 
(>98%) manufactured by 3M. The ESRs were cut into spe- 
cific shapes to allow the coupling of the scintillation light 
from the side of the scintillators. The scintillators were fixed 
to an aluminum frame and interleaved with filters, as shown 
in Fig. 8. 

The scintillation light output is bright when the scintilla- 
tors are irradiated by laser-induced brilliant X-ray pulses, for 
example, >5 x 10” scintillation photons (103 MeV x 54000 
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Table 3. Configuration of the online FSS prototype 


Layer Filter i Sensor Tens 
(mm) (mm) 
1 Al 0.05  Si-PIN 0.3 
2 Al 0.59 Si-PIN 0.3 
3 Al 0 Si-PIN 0.3 
4 Cu 0.11 Si-PIN 0.3 
5 Al 0 Si-PIN 0.3 
6 Al 5.5 Si-PIN 0.3 
7 Cu 45.55 GAGG 2.6 
8 Cu 2.63 GAGG 4.85 
9 Cu 1.3 GAGG 8.67 
10 Cu 11.59 GAGG 8.58 
11 Cu 2.52 GAGG 23.26 
12 Cu 2.06 GAGG 14.26 
13 Cu 12.94 GAGG 9.49 
14 Cu 3.16 GAGG 31.23 
15 Cu 4.91 GAGG 19.58 
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Fig. 8. (Color online) Filters and scintillators used in the prototype 


photons/MeV) per channel per pulse. A non-multiplying pho- 
ton detector is suitable for collecting scintillation light. To 
improve the light-collection efficiency and coupling stability, 
a PIN diode with a large sensitive area (1 cm x 1 cm) was 
used as the photon detector, and the PIN diodes were directly 
coupled with the scintillators using optical grease. As there 
are six PIN diode channels used as X-ray sensors, all PIN 
diodes are mounted on one PIN diode carrier board vertically 
or horizontally, as shown in Fig. 9. 
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Fig. 9. (Color online) PIN diodes and their carrier board 


The sensors, filters, and PIN diodes were placed in an X- 
ray- and electromagnetic-shielding enclosure along with the 
collimator and laser sight, making up the filter and sensor 
unit, as shown in Fig. 10. Bias voltages and electric signals 
were applied or elicited using coaxial cables. 

Using the calibration method described in Section IC, the 
LY x LCE x QE values for the nine GAGG scintillation 
channels were calibrated experimentally, as shown in Fig. 11. 
Error bars are not visible in this figure because of the low peak 
fitting error. 
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Fig. 10. (Color online) Filter and sensor unit 
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Fig. 11. Experimentally calibrated LY x LCE x QE for the nine 
GAGG scintillation channels of the prototype 


B. Electronics unit 


The PIN diode signal is fed to a TIA via alternating cur- 
rent (AC) coupling, and a standard high-speed amplifier with 
a 1.6-GHz gain bandwidth product is adopted as the TIA am- 
plifier. A 2000-ohm resistor connects the cathode of the PIN 
diode to the ground, serving as the direct-current (DC) path, 
which can limit the current of the PIN diode and increase the 
system reliability. The TIA is followed by the main amplifier, 
which adjusts the voltage amplitude to match the dynamic 
range of the DAQ system. In addition to the 16-channel TIAs 
and main amplifiers, the FEE board includes a trigger circuit 
to generate an inner trigger signal for the DAQ system. 

The FEE circuit gains are determined to ensure that the 
signal amplitude is within the dynamic range of each circuit 
stage. First, the energy deposition per pulse in the PIN diodes 
and scintillators was calculated using GEANT4. The spec- 
trum of the X-ray source was modeled using Eq. (9) with a 
photon fluence of 1 x 10!! photons/sr/pulse, and is located 
1 m away from the collimator of the online FSS. Therefore, 
the number of X-ray photons passing through the collimating 
aperture is 2.8 x 10° per pulse. 

The temporal waveforms were calculated using the tran- 
sient response model, which includes the GAGG scintillation 
light, PIN diode output current, and output voltage of each 
circuit stage. The PIN diode was treated as a second-order 
Butterworth low-pass filter (LPF) with a 40-MHz -3 dB band- 
width [58]. The transient response of the GAGG scintillator 
can be expressed as a single exponential decay signal with 
a decay time of 100 ns, and the TIA can also be treated as 
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a second-order Butterworth LPF with a -3 dB bandwidth of 
/GBP/(2xRrCp) [53], where GBP is the 1.6-GHz gain 
bandwidth product, Rp is the feedback resistance, and Cp 
is the input capacitance of the TIA. Cp is mainly attributed 
to the terminal capacitance of the PIN diodes and the stray 
capacitance of the cables, and is in the order of ~100 pF. 
The main amplifier was designed as an inverting second-order 
voltage-controlled voltage source LPF with a 10-MHz -3 dB 
bandwidth. The simulated FEE output waveform obtained by 
the transient response model is plotted along with the exper- 
imental waveform measured in the laser-plasma experiment, 
as shown in Fig. 12. The simulated waveform is in good 
agreement with the experimental waveform, e.g., the relative 
difference in the peak amplitudes is less than 2% when the in- 
tegrals from 130 to 700 ns are equal (with the same energy de- 
positions), thus validating the transient response model built 
in this study. 


— Experimental waveform 
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Fig. 12. (Color online) Simulated and experimental waveforms of 


the FEE output 


As the voltage output swing of the TIA is 0-3 V and the 
input dynamic range of the DAQ system is 1 V, the main am- 
plifier gains are determined to fully utilize the dynamic range 
of the DAQ system, which is 0.32. Then, using the transient 
response model, the TIA gains are determined to ensure that 
the peak amplitudes of the 15-channel TIA outputs are in the 
order of 500 mV. The results are presented in Table 4, with 
the TIA and main amplifier gains listed on the left and right 
sides of the multiplication signs, respectively. The main pur- 
pose of this section is to describe the hardware design. As the 
X-ray fluence can vary significantly in different laser-plasma 
experiments, the FEE is designed as a plug-in module so that 
the gains can be determined using the described methods and 
can be easily changed. 


Table 4. FEE gains 
1-6 7-12 13-15 
1000.32 2000.32 1000x0.32 


Channels 
p in ohm 


To preserve the temporal waveform information and enable 
flexible digital signal processing, a commercial 16-channel 
high-speed analog-to-digital conversion module based on 
DRS4 chips was adopted as the DAQ system. The DAQ sys- 
tem provides a 1-GHz sample rate, a 1024-sample acquisition 
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window (equivalent to 1.024 us at 1 GHz), and a 30-MB/s 
transfer rate (using the USB 2.0 protocol). Given that the 
GAGG decay time is 100 ns, a 1.024-us acquisition window 
is sufficient to capture the complete waveform. The online 
FSS system requires a high repetition rate of several hundred 
kHz. However, the DAQ data transfer rate imposes a limita- 
tion of approximately 1 kHz, calculated as 30 MB/s divided 
by "16 channels x 12 bits/sample x 1024 samples/chan- 
nel/pulse," resulting in 1220 pulses/s. In future campaigns 
in which higher repetition rates may be available, upgrading 
the data transfer protocol would allow for a faster online FSS 
system. 


IV. BREMSSTRAHLUNG EXPERIMENT 


To test the functionality of the online FSS design, the pro- 
totype was implemented in a laser-driven bremsstrahlung ex- 
periment to facilitate the diagnosis of the X/y-ray spectrum. 
Laser-driven bremsstrahlung offers several advantages, in- 
cluding high energy, micro-spot size, and a compact system 
scale, thereby leading to a significant improvement in the 
resolution of high-energy computed tomography (CT) [61]. 
These bremsstrahlung sources typically operate at high rep- 
etition rates. Real-time monitoring of the spectrum plays a 
vital role in evaluating the stability of the radiation source 
and aids in CT image reconstruction. 

The experiment was conducted utilizing a high-repetition 
laser system with a 100-TW power, as shown in Fig. 13. 
During the experiment, a laser pulse was focused onto the 
gas chamber to generate high-energy low-emittance electron 
beams via a laser wakefield acceleration (LWFA) mechanism. 
The laser beam was focused onto a 25-um spot using an off- 
axis parabolic mirror with a 2-m focal length, which captured 
approximately 58% of the laser energy. By reducing the pulse 
duration to as short as 24 fs, a peak intensity of approximately 
4.4 x 10!8 W/cm? was achieved at the target point with an 
energy input of 2 J. A 4-mm-long gas chamber served as 
the target for generating the electron beams. Pure nitrogen 
was used to enhance the plasma density and increase the total 
charge via ionized injection. The molecular density reached 
approximately 1.4 x 1019 cm~? at a pressure of 1000 kPa. 
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Fig. 13. (Color online) Bremsstrahlung experimental setup 


First, the LWFA electron beam properties were assessed 
using an electron spectrometer. The electron beam traversed 
a 300-um-thick beryllium vacuum window, with the elec- 
tron spectrometer positioned downstream of the vacuum win- 
dow. The measured electron spectra are shown in Fig. 14. 
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The LWFA electron beam exhibited a divergence of approx- 
imately 7 mrad and charge of approximately 1000 pC. The 
energy of the electron beam ranged from 5 to 150 MeV. 


Spectrum (pC/MeV) 
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Fig. 14. Spectrum of the LWFA electron beam 


Subsequently, the electron spectrometer was dismounted, : 


and a 3-mm tungsten converter was positioned 25 cm from 


the gas target to generate intense bremsstrahlung radiation. ih 


To deflect the transmission electrons, a 10-cm magnet with 
a magnetic field strength of 0.8 Tesla was placed behind the 
converter. The online FSS, situated at a distance of 1.2 m 


. . 69 
from the tungsten converter, was responsible for detecting E 


the on-axis bremsstrahlung spectrum. The waveform mea- 
surements from a single shot are presented in Fig. 15. For 


r 69. 
the GAGG channels, the waveforms were integrated over a 


time window of 100—700 ns, while for the PIN diode chan- 


nels, integration was performed for 100-250 ns. These inte- - 
grated waveforms constituted the experimental channel data j 
and are denoted by D. There are also some pulse-shaping i 


methods that can be used to calculate D if the X-ray sources 
have a very high repetition rate and the waveforms accumu- 
late [62, 63]. The expected energy deposition (Egep) was 
obtained via Geant4 simulation and the energy-signal coef- 
ficients (ESC) were derived and computed using Eqs. (6) 
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expectation-maximization method with 100 iteration steps. 
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Fig. 15. (Color online) Measured waveforms 
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The bremsstrahlung was generated at a repetition rate of 719 
6st 0.1 Hz, and the real-time responses of the online FSS were 720 
es2 recorded. From the collected waveforms, a subset of five 72 
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Fig. 16. (Color online) Spectra unfolded in the preliminary experi- 
ment 


shots exhibiting suitable pulse amplitudes was selected for 
subsequent unfolding to obtain the bremsstrahlung spectra. 
The average bremsstrahlung spectrum and its standard devi- 
ation for the five shots are depicted by the blue line in Fig. 
16. In addition, the bremsstrahlung spectrum was simulated 
using GEANT4. In the simulation, 5 x 10° electrons were in- 


eso jected into a 3-mm tungsten converter with the electron spec- 


tral shape shown in Fig. 14. The on-axis bremsstrahlung 
within the acceptance angle of the online FSS (approximately 
10~° rad) was recorded. The simulated bremsstrahlung spec- 
trum is depicted by the red line in Fig. 16. The experimental 
spectrum exhibits good agreement with the simulated spec- 
trum, particularly in the energy range of 1—80 MeV. The in- 
cident electrons in the simulation possess a total charge of 
800 pC, which is consistent with the charge measured us- 
ing the electron spectrometer. It is worth noting that the 
bremsstrahlung spectrum from a thick converter is challeng- 
ing to model using analytical formulas owing to the signif- 
icant bremsstrahlung absorption by the converter itself [64]. 
Consequently, the parameter-fitting unfolding techniques em- 
ployed in previous studies [28, 29, 40] are not applicable to 
this particular bremsstrahlung. However, by enhancing the 


racy, the expectation-maximization method can be employed 
to unfold the spectrum with higher precision. The recon- 
structed photon flux reduction caused by self-absorption be- 
low 1 MeV matched well and was consistent with the sim- 
ulated spectrum. Nevertheless, a notable discrepancy exists 
between the simulation and experimental results at approx- 
imately 1 MeV. This discrepancy may be attributed to the 
presence of scattered background in the experiment, which 
was not considered in the Geant4 simulation. 


V. CONCLUSION 


In this paper, we propose an optimization method that uti- 
lizes a genetic algorithm to address the ill-posed nature of the 
FSS system. The condition number of the response matrix 
(RM) is employed as the optimization criterion. To enhance 
the unfolding accuracy across energy bands ranging from tens 
of keV to tens of MeV, we employed high-density GAGG 
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scintillators and silicon PIN diodes. Monte Carlo simula- 
tion results demonstrated that the optimized FSS configura- 
tion yielded significantly improved accuracy when measuring 
typical spectra within this energy range. 

The optimized online FSS featuring PIN diodes as scintil- 
lation photon detectors was successfully implemented. This 
choice enhances the light-collection efficiency and ensures 
stable coupling, further improving the unfolding accuracy. In 
addition, the use of PIN diodes facilitates the quantitative cal- 
ibration of light yields and light collection efficiencies. The 
electrical signals from the PIN diodes were amplified and dig- 
itized using a custom front-end electronics (FEE) module and 
a commercial data acquisition (DAQ) card. The current on- 
line FSS prototype operates at a repetition rate of 1 kHz with 
the potential for further enhancement by increasing the data 
transfer rate of the DAQ card, e.g., using an 80-MB/s opti- 
cal link interface or a USB 3.0 interface. By eliminating the 
need for an image plate (IP) scanner or a high-speed CCD/C- 
MOS, the proposed online FSS system offers the advantages 
of cost-effectiveness and compactness. 

The functionality of the online FSS prototype was vali- 
dated through a laser-driven bremsstrahlung experiment. The 


[1] E. Esarey, C. B. Schroeder, and W. P. Leemans, “Physics of 
laser-driven plasma-based electron accelerators,’ Rev. Mod. 
Phys., 81, 1229(2009). DOI: 10.1103/RevModPhys.81.1229 
E. Esarey, B. A. Shadwick, P. Catravas et at., “Synchrotron 
radiation from electron beams in plasma-focusing chan- 
nels,” Phys. Rev. E, 65, 056505(2002). DOI: 10.1103/Phys- 
RevE.65.056505 

K. Ta Phuoc, S. Corde, C. Thaury et al., “All-optical Comp- 
ton gamma-ray source,” Nat. Photonics, 6, 308(2012). DOI: 
10.1038/nphoton.2012.82 

G. Sarri, D. Corvan, W. Schumaker et al., “Ultrahigh Brilliance 
Multi-MeV gamma-Ray Beams from Nonlinear Relativistic 
Thomson Scattering,’ Phys. Rev. Lett., 113, 224801(2014). 
DOI: 10.1103/PhysRevLett.113.224801 

W. Yan, C. Fruhling, G. Golovin et al., “High-order multipho- 
ton Thomson scattering,” Nat. Photonics, 11, 514(2017). DOI: 
10.1038/nphoton.2017.100 

S. Cipiccia, S. M. Wiggins, R. P. Shanks et al., “A tune- 
able ultra-compact high-power, ultra-short pulsed, bright 
gamma-ray source based on bremsstrahlung radiation from 
laser-plasma accelerated electrons,’ J. Appl. Phys., 111, 
063302(2012). DOI: 10.1063/1.3693537 

S. Corde, K. Ta Phuoc, G. Lambert et al., “Femtosecond x rays 
from laser-plasma accelerators,” Rev. Mod. Phys., 85, 1(2013). 
DOI: 10.1103/RevModPhys.85.1 

F. Albert and A. G. R. Thomas, “Applications of laser 
wakefield accelerator-based light sources,’ Plasma Phys. 
Controlled Fusion, 58, 103001(2016). DOI: 10.1088/0741- 
3335/58/10/103001 

B. Guo, X. Zhang, J. Zhang et al., “High-resolution phase- 
contrast imaging of biological specimens using a stable beta- 
tron X-ray source in the multiple-exposure mode,” Sci. Rep., 9, 
7796(2019). DOI: 10.1038/s41598-019-42834-2 

C. P. Jones, C. M. Brenner, C. A. Stitt et al., “Evaluating 
laser-driven Bremsstrahlung radiation sources for imaging and 


[2] 


[3] 


[4 


= 


[5] 


[6] 


[7 


em 


[8] 


[9] 


[10] 


744 
745 
746 
747 
748 
749 
750 
75 


752 
753 
754 
755 
756 
757 


758 


759 


760 


76 


797 
798 
799 
800 
801 
802 
803 
804 
805 
806 
807 
808 
809 
810 
811 
812 
813 
814 
815 
816 
817 
818 
819 
820 
821 
822 
823 
824 
825 
826 
827 
828 
829 
830 
831 


11 


bremsstrahlung spectrum was recorded in real time and ex- 
hibited good agreement with the simulated bremsstrahlung 
spectrum. In future work, we will conduct additional sys- 
tem test experiments and enhance the electromagnetic com- 
patibility to develop a low-noise system with high unfold- 
ing accuracy. Additionally, we will investigate methodolo- 
gies for identifying improved FSS structures by leveraging 
supercomputers for Monte Carlo RM calculations and opti- 
mizing parameters such as the number of layers. These in- 
vestigations will be conducted with the aim of satisfying the 
accuracy requirements associated with diverse experimental 
scenarios. As the demand for high repetition rates and accu- 
racy increases, this optimized, integrated, compact, and cost- 
effective online FSS technique is promising for advancing 
laser-plasma research and application. 
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